Data assimilation and pollution forecasting in Burgers' equation with model error function

This article presents a correction method for a better resolution of the problem of estimating and predicting pollution, governed by Burgers' equations. The originality of the method consists in the introduction of an error function into the system's equations of state to model uncertainty in the model. The initial conditions and diffusion coefficients, present in the equations for pollution and concentration, and also those in the model error equations, are estimated by solving a data assimilation problem. The efficiency of the correction method is compared with that produced by the traditional method without introduction of an error function.Three test cases are presented in this study in order to compare the performances of the proposed methods. In the first two tests, the reference is the analytical solution and the last test is formulated as part of the “twin experiment”.The numerical results obtained confirm the important role of the model error equation for improving the prediction capability of the system, in terms of both accuracy and speed of convergence.     

Validation of a HPLC/MS method for simultaneous quantification of clonidine,morphine and its metabolites in human plasma

A high‐performance liquid chromatography–tandem mass spectrometry method was developed and validated for the simultaneous quantification of morphine, morphine's major metabolites morphine‐3‐glucuronide and morphine‐6‐glucuronide, and clonidine, to support the pharmacokinetic analysis of an ongoing double‐blinded randomized clinical trial that compares the use of morphine and clonidine in infants diagnosed with neonatal abstinence syndrome. Plasma samples were processed by solid‐phase extraction and separated on an Inertsil ODS‐3 (4 μm) column using an 0.1% formic acid in water–0.1% formic acid in methanol gradient. Detection of the analytes was conducted in the positive multiple reaction monitoring mode. The range of quantitation was 1–1000 ng/mL for morphine, morphine‐3‐glucuronide and morphine‐6‐glucuronide, and 0.25–100 ng/mL for clonidine. Intra‐day and inter‐day accuracy and precision were ≤15% for all analytes across the quantitation range. Extraction recovery rates were ≥94% for morphine, ≥90% for M3G, ≥87% for M6G and ≥ 79% for clonidine. Matrix effect ranged from 85–94% for clonidine to 101–106% for M3G. The method fulfilled all predetermined acceptance criteria and required only 100 μL of starting plasma volume. Furthermore, it was successfully applied to 30 clinical trial plasma samples.     

基于有限体积虚拟区域方法的两相流动直接模拟

为提高计算效率，提出有限体积法离散下的虚拟区域颗粒两相流动直接模拟方法.在控制方程中加入相应的虚拟区域源项，保证了颗粒内部的刚体运动特性.该源项中含有颗粒信息部分及流体信息部分.在每次迭代后，对源项中的流体信息部分进行更新，从而更好地保证颗粒内速度的刚体分布.计算静止颗粒圆柱绕流及单个颗粒的沉降过程，验证了算法的准确性.     

自由基浓度测定方法的再验证及其在煤化学中应用

对前人建立的标准曲线法测煤中自由基浓度进行优化,以DPPH标准样品和基准样品的二次积分面积比值为新参数,结果显示新参数标准曲线法的实测值与理论值相对误差都在5%以内;重复性、复现性实验的相对标准偏差都小于3%。将新参数标准曲线法用于分析不同煤化程度煤和新疆黑山煤(HS)沥青质的自由基浓度,发现随着煤化程度增加,其煤中自由基浓度逐渐增大,从低阶褐煤的8.531×10~(17)/g上升到高阶无烟煤3.37899×10~(19)/g;而在HS煤液化过程中,随着加氢液化温度的升高,其沥青质自由基浓度逐渐下降,从290℃的1.5793×10~(18)/g降到450℃的7.410×10~(17)/g,沥青质自由基浓度变化趋势与其产率变化趋势相一致。     

A Theoretical study on the charge and discharge states of Na-ion battery cathode material,Na1+xFePO4F

Na₂FePO₄F is a promising cathode material for a Na-ion battery because of its high electronic capacity and good cycle performance. In this work, first principle calculations combined with cluster expansion and the Monte Carlo method have been applied to analyze the charge and discharge processes of Na₂FePO₄F by examining the voltage curve and the phase diagram. As a result of the density functional theory calculation and experimental verification with structural analysis, we found that the most stable structure of Na_1.5FePO₄F has the P2₁/b11 space group, which has not been reported to date. The estimated voltage curve has two clear plateaus caused by the two-phase structure composed of P2₁/b11 Na_1.5FePO₄F and Pbcn Na₂FePO₄F or Na₁FePO₄F and separated along the c-axis direction. The phase diagram shows the stability of the phase-separated structure. Considering that Na₂FePO₄F has diffusion paths in the a- and c-axis directions, Na₂FePO₄F has both innerphase and interphase diffusion paths. We suggest that the stable two-phase structure and the diffusion paths to both the innerphase and interphases are a key for the very clear plateau. We challenge to simulate a nonequilibrium state at high rate discharge with high temperature by introducing a coordinate-dependent chemical potential. The simulation shows agreement with the experimental discharge curve on the disappearance of the two plateaus. © 2018 Wiley Periodicals, Inc.     

Kinetics of carbon steel dissolution in ammonium chloride solution containing sodium thiosulfate

The dissolution behavior of carbon steel in ammonium chloride (NH₄Cl) solution containing sodium thiosulfate (Na₂S₂O₃) of various concentrations (0.01 and 0.1 M) was investigated using electrochemical impedance spectroscopy (EIS) and other nonelectrochemical techniques. The weight loss and polarization measurements indicate a significant increase in the NH₄Cl corrosion rate of carbon steel on addition of Na₂S₂O_3. The EIS measurements exhibited two capacitive loops at multiple direct current (dc) potentials for both the concentrations. Electrical equivalent circuit (EEC) and reaction mechanism analysis (RMA) were employed to analyze the impedance data. A four-step mechanism with two intermediate adsorbate species of same charge was proposed to explain the dissolution behavior of carbon steel in the given system. The surface coverage values enumerated that the surface was entirely covered with adsorbed species unlike in the pure NH₄Cl system. Charge transfer resistance and polarization resistance values estimated from RMA parameters indicate the increase in a dissolution rate with dc potential. The surface morphology was inspected via field emission scanning electron microscopy, and the corrosion products including surface state of carbon steel electrode were analyzed using energy dispersive X-ray spectroscopy and X-ray photoelectron spectroscopy.     

The improvement of numerical modeling in the solution of incompressible viscous flow problems using finite element method based on spherical Hankel shape functions

In this paper, the finite element method with new spherical Hankel shape functions is developed for simulating 2‐dimensional incompressible viscous fluid problems. In order to approximate the hydrodynamic variables, the finite element method based on new shape functions is reformulated. The governing equations are the Navier‐Stokes equations solved by the finite element method with the classic Lagrange and spherical Hankel shape functions. The new shape functions are derived using the first and second kinds of Bessel functions. In addition, these functions have properties such as piecewise continuity. For the enrichment of Hankel radial basis functions, polynomial terms are added to the functional expansion that only employs spherical Hankel radial basis functions in the approximation. In addition, the participation of spherical Bessel function fields has enhanced the robustness and efficiency of the interpolation. To demonstrate the efficiency and accuracy of these shape functions, 4 benchmark tests in fluid mechanics are considered. Then, the present model results are compared with the classic finite element results and available analytical and numerical solutions. The results show that the proposed method, even with less number of elements, is more accurate than the classic finite element method.     

Local discontinuous Galerkin method for a nonlocal viscous conservation laws

The purpose of this article is to investigate high‐order numerical approximations of scalar conservation laws with nonlocal viscous term. The viscous term is given in the form of convolution in space variable. With the help of the characteristic of viscous term, we design a semidiscrete local discontinuous Galerkin (LDG) method to solve the nonlocal model. We prove stability and convergence of semidiscrete LDG method in L² norm. The theoretical analysis reveals that the present numerical scheme is stable with optimal convergence order for the linear case, and it is stable with sub‐optimal convergence order for nonlinear case. To demonstrate the validity and accuracy of our scheme, we test the Burgers equation with two typical nonlocal fractional viscous terms. The numerical results show the convergence order accuracy in space for both linear and nonlinear cases. Some numerical simulations are provided to show the robustness and effectiveness of the present numerical scheme.     

基于PI法的门源MS6.4地震前震中附近地震热点图像异常变化研究

基于1980年以来青海省和甘肃省区域地震台网的定位地震目录，分析了研究区最小完备震级，选取空间网格尺度为0.2°×0.2°、地震活动异常学习时间段和预测时间段均为3 a，应用图像信息（PI）法，分析了门源M_S6.4地震发生前震中附近地震热点图像的异常变化过程。结果表明，在2013年1月1日至2016年1月1日的预测时间窗内，冷龙岭断裂和祁连山北缘断裂附近存在明显的地震热点，且地震热点图像颜色偏深，发震概率较高，向前滑动至2014年1月1日至2017年1月1日，2015年1月1日至2018年1月1日和2016年1月1日至2019年1月1日3个预测时间窗，此两断裂附近仍存在地震热点，地震热点图像颜色逐渐由深变浅，而门源M_S6.4地震就落在PI法得到的地震热点内。门源M_S6.4地震前存在明显的地震热点，且地震热点相对集中，地震热点图像颜色由深变浅，震后消失；在祁连山北缘断裂附近，地震热点持续存在，分布较集中，地震热点图像颜色也呈由深变浅的过程。未来需关注祁连山北缘断裂，附近有发生强震的可能。